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CALCULATION  OF  MINIMUM  ENTRY  HEAT  TRANSFER  SHAPE 
OF  A  SPACE  VEHICLE 

Zhou  Qi  cheng 

ABSTRACT 

This  paper  dealt  with  the  study  of  the  minimum  total 
heat  transfer  aerodynamic  shape  problem  for  a  re-entry 
spacecraft.  In  this  paper,  the  well  known  boundary  layer 
heat  transfer  relation  was  used  first  to  derive  the 
expressions  of  the  total  heat  flow  rate  on  the  surface  of 
the  spacecraft  and  the  total  heat  added  during  the  entire 
re-entry  process.  Then,  with  the  aid  of  some  numerical 
results,  the  flow  field  and  the  Prandtl- Meyer  expansion 
theory,  the  minimum  re-entry  heat  transfer  shape  under 
specified  fineness  ratio  and  total  vehicle  weight  conditions 
could  be  obtained  using  a  variational  method.  Finally, 
the  effect  of  the  drag  coefficient  on  the  minimum  re-entry 
heat  transfer  shape  was  analyzed.  The  calculated  results 
indicated  that:  For  a  given  total  vehicle  weight  and  spe¬ 
cified  fineness  ratio,  the  total  heat  transfer  decreased 
with  increasing  drag  coefficient.  The  minimum  re-entry  heat 
transfer  shape  is  a  cylinder  with  a  flat  nose. 

I .  INTRODUCTION 

The  rational  section  of  the  aerodynamic  shape  of  the  re-entry 
vehicle  to  reduce  the  heating  due  to  surface  convection  to  the  maximum 
extent  possible  is  very  meaningful  to  the  design  of  either  retrievable 
satellites  or  missile  warheads.  It  not  only  can  effectively  reduce 
the  weight  of  the  re-entry  heat  protection  layer  to  improve  the 
effective  payload  of  the  vehicle,  but  also  is  helpful  in  improving 
the  accuracy  of  the  calculation  of  the  aerodynamic  coefficients  to 
improve  the  target  hitting  rate.  Therefore,  in  recent  years,  many 
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researchers  and  engineering  technical  personnel  are  concerned  about 
this  problem.  This  paper  involved  the  study  of  the  minimum  heat 
transfer  aerodynamic  shape  problem  of  a  re-entry  vehicle  on  the  basis 
of  the  design  of  a  retrievable  satellite.  In  order  to  simplify  the 
problem,  we  assumed  that  during  the  re-entry  flight  process  the 
geometrical  shape  of  the  vehicle  remained  unchanged.  The  geometrical 
slope  of  the  surface  of  the  object  is  greater  or  equal  to  zero. 

References  [1,2]  studied  the  minimum  heat  flow  rate  shape  with 
a  specified  fineness  ratio  when  the  flight  Mach  number  remained  un¬ 
changed  under  the  assumption  of  Newtonian  pressure.  Reference  [3] 
extended  the  problem  to  the  entire  re-entry  flight  and  discovered 

that  for  a  given  trajectory  coefficient/-^  \and  a  specified  fineness 

V  C„A  / 

ratio,  the  minimum  heat  transfer  shape  is  the  minimum  heat  flow 

rate  shape.  Reference  [4]  applied  the  numerical  values  of  inviscid 
flow  to  the  problem  discussed  in  [3].  The  accuracy  of  the  calculation 
of  the  minimum  heat  transfer  shape  was  greatly  improved.  The  minimum 
heat  transfer  shape  (which  is  also  the  minimum  heat  flow  rate  shape) 
under  laminar  and  turbulent  flow  conditions  for  a  given  trajectory 
coefficient  obtained  in  that  paper  was  a  rotating  body  with  a  flat 
nose.  The  angles  of  inclination  at  the  tip  and  the  surface  were 
about  20°  and  17°,  respectively.  In  [5],  the  actual  requirement  of 
an  effective  volume  In  the  vehicle  was  taken  Into  account  during  the 
engineering  design.  Therefore,  the  work  carried  out  in  [3,4]  was  fur¬ 
ther  extended.  The  minimum  heat  transfer  shape  of  a  re-entry  flight 
vehicle  in  laminar  and  turbulent  flow  condition  with  a  given  trajec¬ 
tory  coefficient,  a  specified  fineness  ratio  and  a  specified  effective 
volume  was  calculated.  Furthermore,  the  analytical  expression  of  the 
minimum  heat  transfer  shape  under  laminar  flow  condition  with  a  speci¬ 
fied  trajectory  coefficient  and  a  specified  effective  volume  (no 
requirement  was  needed  for  the  fineness  ratio)  was  obtained.  The 
above  references  did  not  involve  the  effect  of  the  variation  of  drag 
coefficient  on  the  minimum  heat  transfer  shape. 

In  the  actual  engineering  design,  people  frequently  wanted  to 
know  the  shape  of  the  minimum  re-entry  heat  transfer  under  a  given 
total  weight  and  a  specified  fineness  ratio.  Because  the  re-entry 
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trajectory  coefficient  is  an  Important  parameter  influencing  the 
re-entry  flight,  the  variation  of  drag  coefficient  will  directly 
vary  the  trajectory  coefficient.  Consequently,  the  total  added  heat 
to  re-entry  will  be  affected.  Therefore,  in  the  optimization  of  the 
aerodynamic  shape,  it  is  necessary  to  analyze  the  effect  of  the  drag 
coefficient  of  the  re-entry  body  on  the  re-entry  heat.  The  study  of 
whether  the  minimum  heat  flow  rate  shape  is  still  the  minimum  re-entry 
heat  transfer  shape  is  very  important.  The  purpose  of  this  paper  is 
to  find  the  minimum  re-entry  heat  transfer  shape  with  specified  fine¬ 
ness  ratio  and  total  vehicle  weight. 

II.  PRESENTATION  OF  THE  PROBLEM 


The  commonly  adopted  LEES  laminar  flow  heat  transfer  theory  [6] 
and  VAGLIO-LAURIN  turbulent  flow  heat  transfer  theory  [7]  are  used 
to  calculate  the  surface  heat  flow  rate  of  the  re-entry  vehicle 


where 


0.0332Pr~,,,g(  1  —  Hm/Ht) 
0 . 0296 1  -  HJH%) 


laminar  flow 
turbulent  flow 


p,  u,  y  and  H  represent  the  density,  tangential  velocity,  viscosity 
coefficient  and  enthalpy  of  the  flow,  respectively 

(p* “ P./p». t,— «,/ y/  tH»  ,  }  y  and  s  represent  the  dimensionless 

abscissa  of  the  vehicle  shape  and  the  arc  length  along  the  surface, 
respectively  ^*)  ,  where  Rg  is  the  radius  of  the  base. 

Pr  is  the  Prandtl  number  and  g  is  the  gravitational  acceleration. 

The  subscripts  e,  o  and  w  represent  the  values  at  the  outer  fringe, 
stationary  point,  and  the  wall  of  the  boundary  layer,  respectively, 
n  »  2  represents  laminar  flow  and  n  =  5  represents  turbulent  flow. 

By  introducing  the  following  symbol  @ 

9-p, *- 1 

and  integrating  the  heat  flow  rate  along  the  surface  of  the  flight 
vehicle,  we  can  obtain  the  total  heat  flow  rate  on  the  surface  of  the 
flight  vehicle  as 


rf. 

UTDfJJZF*—  LJ.^'J 


It  is  commonly  known  that  the  stationary  point  heat  flow 
rate  is  related  to  the  radius  of  curvature  of  the  object  at  the 
stationary  point.  The  magnitude  of  the  heat  flow  rate  increases  or 
decreases  with  increasing  or  decreasing  radius  of  curvature.  There¬ 
fore,  it  is  reasonable  to  believe  that  the  minimum  heat  transfer 
shape  is  composed  of  a  flat  nose  and  a  rear  body  part  [1,4].  Hence, 
by  using  the  isoenthalpy  relation  and  the  viscosity  exponent  law, 
we  can  express  the  integral  in  the  above  equation  as 


where 


■ 

*'(*-)-[  (y_*)Jl#r  +  l] 

a_  l+(v-l)fi> 
v 


Y  is  the  specific  heat  of  air,  to  is  the  viscosity  exponent,  y^  is 

the  height  of  the  flat  nose,  is  the  flight  Mach  number  and  P  is 

°°  P  Ax 

the  pressure  coefficient  which  is  defined  as  is  the 

r%  ay  y 

dimensionless  ordinate  of  the  flight  vehicle  shape  and  I. 

Al 

is  the  integral  on  the  flat  nose. 


If  we  assume  that  the  sonic  point  is  located  at  the  intersection  of 
the  flat  nose  and  the  rear  body,  then  it  is  possible  to  complete  this 
integration  by  using  the  simplified  integral  relation  [8]  or  an  inviscid 
flow  numerical  value.  Therefore,  it  is  a  known  constant.  We  can 
rewrite  it  as 


»_/  v— 1  r  ,+-r 

*'  \y+r)  /•  * 


where  is  a  known  constant  and  its  value  differs  for  laminar  flow 
and  turbulent  flow. 


Based  on  the  differential  correlation  between  the  total  heat 
flow  rate  Q  and  the  total  re-entry  heat  Q 


■  r  f*i  dt 

0-j  0-^* 


and  the  CHAPMAN  [9]  trajectory  approximation  which  is  applicable  to 
trajectories  with  larger  variation  of  the  angle  of  inclination  6 

dt 

we  can  express  the  total  re-entry  flight  heat  as 

Q— *pr[J>  F~ 

where  Vot  is  the  flight  speed,  h  is  the  flight  altitude,  t  is  the 
flight  time  and  h£  is  the  upper  fringe  height  of  the  dense  atmosphere. 

For  an  Earth  orbit  re-entry  flight  vehicle,  aerodynamic  heating 
primarily  occurs  during  the  hypersonic  flight  stage  (M-»l)  .  in  an 

approximation  analysis,  it  can  be  considered  that  the  pressure  coeffi¬ 
cient  of  the  rear  body  is  not  related  to  the  Mach  number  of  the 
incoming  flow  (such  as  the  Newtonian  theory  or  the  P-M  expansion 
theory).  It  is  only  a  function  of  the  geometric  position.  Therefore, 
the  integration  with  respect  to  <J>  in  the  above  formula  can  be  ex¬ 
tracted  from  the  integration  with  respect  to  h.  Hence,  it  is  possible 
to  use  the  relationship  before  and  after  the  normal  shock  wave,  the 
isoenthalplc  correlation  and  the  CHAPMAN  approximation  trajectory  to 
finally  simplify  the  total  re-entry  heat  formula  under  the  condition 
thatA/„»l  as 


where 


’-Mot]1’- 

,r  u  1  *<*->  d 
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Figure  1.  The  variation  of  B  with  the  re-entry  angle 
and  the  trajectory  coefficient 
1 — laminar  flow;  2 — turbulent  flow 

8  is  the  reciprocal  of  the  density  of  the  atmosphere,  and  R  is  the 
average  radius  of  the  Earth.  From  [91,  we  know  that  the  dimension¬ 
less  CHAPMAN  variables  Z  and  dZ/dU  are  known  functions  of  the  dimen¬ 
sionless  velocity  U  which  is  only  related  to  the  re-entry  angle  0„. 

hj 

The  functional  relationship  between  M^  and  U  is  only  related  to 
M/CpA.  Actual  computation  showed  that  B  is  a  function  of  the  re-entry 
angle.  The  effect  of  the  trajectory  coefficient  can  be  completely 
neglected.  The  details  are  shown  in  Figure  1.  A  is  the  reference 
area  of  the  drag  coefficient  CD,  M  is  the  mass  of  the  flight  vehicle 
and  A  is  a  constant.  Its  expression  is 

m  _  i  — !_  t  — —  l  — -i— +-2L. 

Am-- ^-2  •  C„mR,  ■  (0R)  a  u  * 

- ,  +  1  1  +  <»  1  tm 

X  0  »■  R‘  '  HI*  v.  ' 

where  Vg  and  are  the  re-entry  velocity  and  the  air  viscosity 
coefficient  at  the  initial  re-entry  altitude. 

T'  a  so-ca'1  ed  minimum  heat  transfer  shape  is  the  geometric  shape 
with  whxcr  tl.  total  re-entry  heat  Q  of  the  flight  vehicle  surface  is 
the  minimum.  If  x  *  x(y)  is  used  to  represent  the  equation  of  the 
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parent  line  of  the  minimum  heat  transfer  shape,  it  can  be  substituted 
into  the  computational  expression  of  Q  derived  above  as  equation  (1) 
and  Q  will  have  the  minimum  value.  Therefore,  solving  for  the  min¬ 
imum  heat  transfer  shape  under  specified  fineness  ratio  and  total 
vehicle  weight  can  be  mathematically  reduced  to  the  solving  of  the 
minimum  curve  of  the  generalized  function  I/CD  inside  the  bracket  in 
(1).  That  is  to  find  the  minimum  curve  of  the  following  mixed  type 
generalized  function  with  moving  boundary  points. 


/•[*(>)]  = 


where 


/v— 1\**»  .  1  +  »- 1  1 1 

(•JT+ij  1  KVi  +J  ,,F(P,±,y)dy 


4  J,  •  pyJy+  *\ytpydx 


The  first  term  of  the  formula  represents  the  contribution  of  the 
flat  nose  to  the  drag.  It  is  possible  to  use  the  simplified  integral 
relation  [8]  to  complete  this  integration.  Or,  it  is  also  possible 
to  compute  it  based  on  the  results  of  the  numerical  values  of  an 
inviscid  flow.  The  second  term  represents  the  contribution  of  the 
rear  body  to  the  drag.  In  order  to  facilitate  the  use  in  the  future, 
we  arc  rewriting  the  above  formula  as 

Co  *■  4/ By [  +  4  Pydy  (  3 ) 

where  JD  is  a  known  constant. 

If  t  is  used  to  represent  the  fineness  ratio(r »tRa/L),  of  the 
flight  vehicle,  then  the  boundary  condition  can  be  expressed  as 

2  (4) 

Because  we  did  not  specify  the  height  of  the  flat  nose  y^,  it 
is  a  quantity  to  be  determined.  In  order  to  have  a  solution  to  the 
variational  problem,  we  must  supplement  the  so-called  crossing  condi¬ 
tion  to  be  satisfied  by  the  minimum  curve  at  x  =  0.  Its  actual  form 
is  related  to  the  treatment  method  of  the  problem. 


Equations  (2),  (4)  and  the  crossing  condition  at  x  =  0  are  the 
original  equations  to  solve  for  the  minimum  re-entry  beat  transfer 


shape  for  a  specified  total  vehicle  weight  and  fineness  ratio. 

III.  SOLUTION  OF  THE  VARIATIONAL  EQUATION 

The  minimum  heat  transfer  shape  for  a  given  total  vehicle  weight 
M  and  specified  fineness  ratio  t  can  be  found  by  solving  the  varia¬ 
tional  problem  in  equations  (2)  and  (4)  directly.  However,  in  order 
to  facilitate  the  analysis  of  the  effect  of  the  variation  of  drag 
coefficient  Cp  on  the  minimum  heat  transfer  shape,  we  chose  to  use 
the  following  method.  It  involved  the  use  of  a  series  of  properly 
chosen  specified  values  of  drag  coefficients  C0<(i— 1,2,3—)  .  Then,  the 

minimum  beat  transfer  shapes  and  the  total  re-entry  beats  correspond 
to  these  specified  values  of  drag  coefficients.  Therefore,  we 
can  understand  the  effect  of  drag  coefficient  on  the  minimum  heat 
transfer  shape  and  the  total  heat  from  this  series  of  calculations. 

In  the  meantime,  it  is  not  difficult  to  find  the  eventually  needed 
minimum  heat  transfer  shape  for  a  given  total  weight  and  a  specified 
fineness  ratio. 


According  to  our  knowledge  in  the  variational  method  [10],  the 
minimum  re-entry  heat  transfer  shape  for  a  specified  drag  coefficient 
is  reduced  to  the  solution  of  the  minimum  curve  for  the  following 


generalized  function 

/**C*Cy)]-(-^-)*'Vjry1  "  1  -X,CD+U,Joy\  +  J  ,F**dy 

y.  i 

where  X1  is  the  Lagrange  coefficient  yet  to  be  determined 

F**-F+iX,Py 

The  boundary  condition  is  the  same  as  before  as  shown  in  equation  (4) 
The  corssing  condition  is 

H1 + irrr>'  +«vw-.-(. )wi  -  o  (6) 

In  order  to  determine  the  Lagrange  coefficient  X^,  it  is  necessary 
to  simultaneously  solve  for  the  constraint  equation  (3)  in  which  the 
drag  coefficient  is  equal  to  a  specified  value. 


Up  until  the  present  moment,  we  only  assumed  that  the  junction 
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Figure  2.  The  minimum  heat  transfer  shape  for  a  given 
trajectory  coefficient  (t  =  1.0) 

2 — laminar  flow;  3 — turbulent  flow;  4 — xaminar  flow;  5 — turbulent 
flow;  6 — this  paper;  7 — reference  [4];  8 — reference  [3] 

between  the  flat  nose  and  the  rear  body  is  a  sonic  point  and  its 
pressure  coefficient  is  not  related  to  the  Mach  number  of  the  incoming 
flow.  References  [1-3]  used  the  Newtonian  theory  to  calculate  the 
pressure  distribution  on  the  rear  body.  However,  the  minimum  heat 
transfer  shape  they  obtained  is  quite  different  from  the  one  obtained 
in  [4]  using  the  inviscid  flow  numerical  solution  as  a  basis.  By 
taking  into  account  that  fact  that  the  variational  method  requires 
that  the  rear  body  pressure  distribution  starting  from  the  sonic 
point  must  be  continuous  and  differentiable,  we  used  the  P-M  expansion 
theory  to  calculate  the  rear  body  pressure  distribution  in  this  paper. 
Although  by  doing  so  it  is  possible  that  the  numerical  value  of  the 
rear  body  pressure  coefficient  might  be  underestimated,  yet  in  actual 
trial  computation,  it  was  discovered  that  the  major  factors  affecting 
the  minimum  heat  transfer  shape  are  the  pressure  distribution  on  the 
flat  nose  and  the  pressure  variation  near  the  junction  between  the 
flat  nose  and  the  rear  body.  By  using  the  P-M  expansion  theory,  the 
pressure  variation  near  the  junction  can  be  better  reflected.  From 
Figure  2,  we  can  see  that  the  minimum  heat  transfer  shape  for  a  given 
trajectory  coefficient  using  the  P-M  expansion  theory  in  this  paper 
agrees  with  that  obtained  based  on  the  inviscid  flow  numerical  solu¬ 
tion  in  [M]  comparatively  well.  As  can  be  expected,  when  the  calcul¬ 
ated  results  converged  to  approach  the  minimum  curve  (later,  we  will 
see  that  the  minimum  heat  transfer  shape  is  a  cylinder  with  a  flat 
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nose),  the  accuracy  of  the  rear  body  pressure  distribution  calculated 
using  the  P-M  expansion  method  is  improved.  In  addition,  we  should 
also  point  out  that  the  treatment  in  this  paper  is  focused  on  the 
study  of  the  variation  trend  of  the  target  function  (i.e.,  the  total 
re-entry  heat)  rather  than  the  numerical  value  of  the  target  function 
itself.  The  computational  error  of  the  rear  body  pressure  distribu¬ 
tion  will  not  cause  the  variation  of  the  trend  of  the  variation  of  the 
total  re-entry  heat.  Therefore,  it  is  believed  in  this  paper  that  it 
is  feasible  to  use  the  P-M  expansion  theory  to  calculate  the  rear  body 
pressure  distribution  as  an  approximate  analysis. 

When  we  use  the  P-M  method  to  calculate  the  pressure  on  a  conical 
body,  the  pressure  coefficient  is  only  related  to  the  local  slope  dx/dy. 
According  to  our  knowledge  in  the  variational  method,  the  Euler  equa¬ 
tion  of  this  type  of  variational  problem  has  a  simple  integral: 

-fr+4Ao'-|f-c 

where  C  is  an  integration  constant  yet  to  be  determined  by  the  bound¬ 
ary  condition.  The  above  equation  can  be  integrated  once  more  into 
the  following  form: 

r"  u| 

- -  ZET - (7) 

n  ....  «-i  ,  JUL 

7^rrA*)y  -4A,  at 


yrrT-/(*)-4A.|f - C 


where 


and  *2  are  the  slopes  of  the  minimum  curve  at  x 
respectively.  They  are  both  unknown  quantities. 


=  0  and  x  =  2/t 


Substituting  the  correlations  of  the  P-M  expansion  theory 
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into  the  above  equation,  it  is  possible  to  find  the  actual  computa- 
tional  expressions  for and  through  not  too  complicated 

derivations.  Therefore,  by  simultaneously  solving  for  equations  (3), 
(4),  (6)  and  (7),  we  can  obtain  the  minimum  re-entry  heat  transfer 
shape  for  a  drag  coefficient  equal  to  a  specific  value. 

IV.  COMPUTATIONAL  RESULTS  AND  DISCUSSION 

We  computed  the  minimum  re-entry  heat  transfer  shapes  for  various 
fineness  ratios  with  drag  coefficients  as  a  series  of  specific  values 
under  the  conditons  that  v»1.2,o-l.o.  and  dy/dx>0  .  The  numerical 
values  of  JK  and  JD  were  obtained  from  the  results  in  [4],  In  the 
following,  we  will  briefly  introduce  the  calculated  results. 

Figure  3  plots  out  the  minimum  re-entry  heat  transfer  shapes 
with  drag  coefficients  equal  to  a  series  of  specific  values.  CD0  is 
the  drag  coefficient  of  the  minimum  heat  flow  rate  shape  (i.e.,  the 
minimum  heat  transfer  shape  studied  in  [4]).  From  this  figure,  we 
can  see  that  the  minimum  re-entry  beat  transfer  shape  is  a  family  of 
flat  headed  rotating  bodies.  With  increasing  drag  coefficient,  the 
height  of  the  flat  nose  is  gradually  higher  and  the  half  vertex  angle 
of  the  rear  body  is  becoming  smaller.  The  rear  body  more  and  more 
approaches  a  cone.  The  laminar  flow  minimum  re-entry  heat  transfer 
shape  and  the  turbulent  flow  minimum  re-entry  heat  transfer  shape  are 
very  close. 

Figure  4  plots  the  minimum  re-entry  beat  transfer  shape  total  beat 

flow  rate  for  various  drag  coefficient  equal  to  different  specific 
0 

values.  Qq  is  the  minimum  heat  transfer  shape  total  heat  flow  rate 
under  invariant  trajectory  coefficient  conditions.  From  the  figure, 
we  can  see  that  when  the  drag  coefficient  is  not  equal  to  C^q,  the 
total  heat  flow  rate  Q  is  always  higher  than  Q  .  It  is  true  for  both 
laminar  flow  and  turbulent  flow.  The  reason  is  very  simple. 
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Figure  3*  The  minimum  re-entry  heat  transfer  shape  when 
the  given  drag  coefficient  is  a  fixed  value 
2 — turbulence;  3 — laminar  flow;  4 — turbulent  flow 


2 — turbulent  flow;  3 — laminar  flow;  4 — turbulent  flow 
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Figure  4.  The  variation  of  the  minimum  heat  transfer 
shape  total  heat  flow  with  drag  coefficient  when  the  drag 
coefficient  is  a  fixed  value 

1 — turbulent  flow;  2 — laminar  flow;  3 — calculated  results  in  the 
paper;  4 — flat  heated  cone  results  from  [4];  5 — turbulent  flow; 

6 — laminar  flow 

It  has  been  pointed  out  in  [4,5]  that  when  the  trajectory  coefficient 
is  unchanged,  the  minimum  re-entry  heat  transfer  shape  is  the  minimum 
beat  flow  rate  shape.  Therefore,  when  Cp  is  different  from  CDQ,  the 
Q  is  naturally  higher  than  Qq. 

Figure  5  plots  the  minimum  re-entry  heat  transfer  shape  total 
heat  corresponding  to  various  drag  coefficients  which  are  equal  to 
different  specific  values.  The  reference  values  used  in  the  figure 
are  similar  to  those  in  Figure  4.  From  this  figure,  we  can  clearly 
see  that  the  re-entry  total  heat  rapidly  decreases  with  increasing 
drag  coefficient.  When  the  drag  coefficient  CD  increased  to  the  value 
equal  to  that  of  a  flat  nosed  cylinder,  the  re-entry  total  beat 
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reached  the  minimum.  This  variation  trend  indicates  that }although 
the  increase  in  drag  coefficient  would  cause  the  surface  total  heat 
flow  rate  to  increase,  yet,  under  a  given  total  re-entry  weight  con¬ 
dition,  the  increase  in  drag  coefficient  would  simultaneously  cause 
the  trajectory  coefficient  to  drop  rapidly.  The  net  effect  is  to 
gradually  reduce  the  total  heat.  From  Figure  5,  it  is  not  difficult 
to  see  that  when  the  re-entry  weight  is  specified,  the  minimum  re-entry 
heat  transfer  shape  is  a  flat  headed  cylinder. 

Figures  4  and  5  also  plotted  the  calculated  results  of  the  flat 
headed  cone  part  using  the  inviscid  flow  numerical  solution  as  the 
basis.  The  variation  trend  agrees  with  the  calculated  results  in  this 
paper.  When  the  drag  coefficient  is  identical,  the  total  heat  flow 
rate  and  the  total  heat  of  the  flat  headed  cone  are  higher  than  the  cal¬ 
culated  results  of  the  minimum  re-entry  heat  transfer  shape  in  this 
paper. 

Figure  6  plots  the  calculated  results  of  a  flattened  cone  and  a 
spherical  cone  which  are  very  close  to  the  minimum  heat  transfer  shape 
[4].  From  the  figure,  we  can  see  that  the  drag  coefficient  of  a  flat 
headed  cone  is  higher  than  that  of  a  spherical  cone  with  the  same 
half  conical  angle.  When  their  half  conical  angles  are  identical, 
although  their  total  heat  flow  rates  (i.e.,  I)  are  not  too  different, 
yet  the  re-entry  total  heat  of  the  flat  nosed  cone  (i.e.,  I/Cp)  is 
much  lower  than  that  of  the  spherical  cone.  In  addition,  we  can  also 
see  from  Figure  6  that  In  the  half  conical  angle  range  shown,  regard¬ 
less  of  the  spherical  cone  or  the  flat  headed  cone,  the  surface  total 
heat  flow  rate  gradually  increases  with  Increasing  drag  coefficient 
(I.e.,  the  re-entry  body  becomes  dull).  The  re-entry  total  heat  gra¬ 
dually  decreases. 

According  to  the  calculated  results  shown  in  Figures  4,  5  and  6 
we  seem  to  be  able  to  draw  the  following  conclusions:  for  a  speci¬ 
fied  total  re-entry  weight,  during  the  selection  of  the  shape  of  the 
re-entry  vehicle,  we  should  choose  a  flat  headed  rotating  body  on  the 
basis  of  satisfying  other  objectives  of  the  whole  body  in  order  to 
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Figure  5.  The  variation  of  minimum  beat  transfer  shape 
total  heat  with  drag  coefficient  when  the  drag  coefficient 
is  a  fixed  value 

1 — calculated  results  in  this  paper;  2 — the  flat  headed  cone  results 
in  [4];  3 — laminar  flow;  ^ — turbulent  flow;  5 — turbulent  flow; 

6 — laminar  flow 
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Figure  6.  Comparison  of  heating  of  a  flat  headed  cone  to 
a  spherical  cone  (t  =  1.0) 

1 — flat  headed  cone;  2 — spherical  cone;  3 — degree 
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Figure  6  (continued) 

1 — turbulent  flow;  2 — laminar  flow;  3 — degree;  4 — laminar  flow; 

5 — turbulent  flow;  6 — degree 

greatly  reduce  the  re-entry  heat  to  the  maximum  extent  possible  and 
to  reduce  the  weight  of  the  re-entry  heat  proof  layer.  Furthermore, 
we  should  attempt  to  use  a  larger  drag  coefficient  to  the  extent 
possible.  The  flat  headed  cone  is  an  aerodynamic  shape  with  some 
future. 


V.  BRIEF  CONCLUSIONS 
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1.  The  minimum  re-entry  beat  transfer  shape  when  the  drag  coeff¬ 
icient  is  equal  to  some  specific  number  is  a  family  of  flat  beaded 
rotating  bodies. 

2.  Along  with  increasing  drag  coefficient,  the  total  beat  flow 
rate  of  the  flat  headed  rotating  body  is  reduced  monotonically .  After 
reaching  the  minimum,  it  monotonically  increases.  The  re-entry  total 
heat  monotonically  decreases  with  increasing  drag  coefficient. 


ft 
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3.  The  minimum  re-entry  heat  transfer  shape  for  a  given  re-entry 
total  weight  is  different  from  the  minimum  beat  transfer  shape  for  a 
given  trajectory  coefficient.  It  Is  a  flat  headed  cylinder. 

4 .  In  order  to  reduce  the  re-entry  heat  to  the  maximum  extent 
possible,  while  satisfying  other  functional  objectives,  we  should 
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choose  a  flat  head  rotating  body.  Furthermore,  a  larger  drag  coeff¬ 
icient  shape  should  be  used  to  the  extent  possible.  The  flat  headed 
cone  is  a  promising  aerodynamic  shape. 
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LIMIT  ANALYSIS  OF  THE  STRUCTURAL  STRENGTH  FOR  THE  COMBUSTION 
CHAMBER  OF  A  LIQUID  ROCKET  ENGINE 


ABSTRACT 

This  paper  is  an  actual  example  of  the  limit  anal¬ 
ysis  method  in  mechanics.  It  provided  computational  for¬ 
mulas  for  the  limiting  load  of  the  major  parts  in  the  com¬ 
bustion  chamber  of  a  liquid  rocket  engine.  These  calculated 
results  have  good  engineering  accuracy  in  approximation 
through  experimental  verification. 


I .  INTRODUCTION 

Due  to  the  "short  time"  and  "single  use"  characteristics  of  the 
liquid  rocket  engine,  it  is  believed  to  be  reasonable  for  part  of  its 
components  to  work  in  a  plastic  state  or  for  the  entire  unit  to 
approach  the  plastic  state.  To  design  a  liquid  rocket  engine  accord¬ 
ing  to  this  viewpoint  will  enable  us  to  fix  the  load  resistance  poten¬ 
tial  of  the  structure  to  the  maximum  degree  possible.  Consequently, 
the  weight  strength  ratio  of  the  engine  can  be  rationally  designed. 

In  the  60's,  the  workers  in  mechanics  in  our  country  con¬ 
ducted  a  lot  of  work  in  the  areas  of  limit  analysis  and  limit  design 
to  solve  the  common  engineering  shell  limit  analysis  method.  Further¬ 
more,  this  method  was  developed  to  establish  its  theoretical  basis  to 
promote  the  flourishing  development  of  this  technology.  This  situation 
makes  it  not  only  necessary,  but  also  possible  to  use  the  limit  anal¬ 
ysis  method  in  a  liquid  rocket  engine. 

The  limit  analysis  method  is  an  important  branch  in  plastic 
mechanics.  Rigorously  speaking,  it  is  an  engineering  approximation 
method.  Actual  practice  proved  that  the  conclusions  obtained  by  using 
this  method  had  the  sufficient  engineering  approximation.  Because  of 
the  ease  of  use  and  the  clear  physical  concept  of  this  method,  we 

This  paper  was  received  on  July  31,  1981. 
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believe  that  the  application  of  limit  analysis  to  the  design  of 
liquid  rocket  engine  strength  should  have  some  constructive  effect 
for  a  long  time  to  come. 

Finally,  it  should  be  pointed  out  that  all  the  components  of  the 
liquid  rocket  engine  are  located  in  a  strong  vibration  environment. 
Dynamic  characteristic  design  and  calibration  are  sometimes  necessary 
for  some  components.  This  paper  did  not  involve  any  discussion  in 
this  aspect. 

II.  THE  LIMIT  ANALYSIS  OF  THE  LOADED  TOP  COVER  OF  THE 
COMBUSTION  CHAMBER 


The  form  of  the  loaded  top  cover  is  dependent  on  the  supporting 
condition  of  the  combustion  chamber.  It  is  also  related  to  the  re¬ 
quired  dimension  of  the  entire  unit.  The  usual  forms  are  the  spheri¬ 
cal  top  cover,  which  is  a  part  of  the  sphere.  It  is  a  more  advantag¬ 
eous  structural  form  to  sustain  a  uniformly  distributed  pressure. 

The  triple  centered  top  cover.  By  using  such  a  structural  form, 
it  is  not  only  possible  to  take  care  of  the  favorable  form  in  main¬ 
taining  the  structure  under  high  working  pressure,  but  also  able  to 
take  the  requirement  of  reducing  the  total  length  dimension  of  the 
entire  unit  into  account. 


The  conical  top  cover.  When  the  thrust  of  the  engine  is  trans¬ 
ferred  from  the  head  of  the  combustion  chamber,  the  conical  cover  top 
is  the  most  favorable  force  exerting  component. 


1.  Spherical  top  cover  with  an  opening  in  the  center  to  sustain 
uniform  pressure.  The  simplified  mechanical  model  of  the  structure 
is  shown  in  Figure  2-1.  The  internal  boundary  is  connected  to  the 
rigid  flange.  The  external  boundary  is  considered  as  the  solid  sup¬ 
port.  When  the  limiting  condition  is  reached,  let  us  assume  that  the 
distortion  velocity  is 


W ~  -  A(  coup,—  co*<p) 


(2-1) 
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As  for  the  direction  of  W,  the  outer  normal  direction  is 
positive . 


Figure  2-1 


Using  the  mobile  method  and  the  jacketed  shell  yield  condition, 
we  can  get ,  when  *•>*,. 

tin’?,  (l—K)(ain*p*—tin,p,)+2Ktin,p,+2(cosp0  —  coBp,)cos<p, 

3  K  i  (  rn*  m .  —  roxm.  Yf  rni  m  —  rt\*  .  X  rin  ^  0 


when  *•<*>> 


where 


(cos??,—  cosf),)(co«9?,  —  co*«5,+  tin  *<p, 
_ _ 2  sin  *p. 


3(cos??i—  cosip, )(  cos  9?,  —  cos??, +  tin  *<p , 


(2-2) 


(2-3) 


(2-4) 


(2-5) 


2  P/?’  sin  *?», 

3 


(2-6) 


The  numerical  results  obtained  from  equations  (2-2)  and  (2-3) 
are  tabulated  in  Table  2-1. 


The  limiting  loading  capability  curves  are  plotted  as  shown  in 
Figure  2-2. . 
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TABLE  2-1.  The  limiting  loading  capability  of  an  opened  spherical 
shell  with  external  boundary  solid  support. 


1 

[  *.-0' 

6»,-10* 

«»,-20* 

6». 

Jf-0.01 

0.02 

o.os 

0.01 

0.02 

0.05 

0.01 

0.02 

0.05 

10* 

2.961 

JO* 

6.171 

2.969 

2.640 

4.150 

3.112 

2.099 

so* 

9.606 

6.429 

9.007 

9.009 

4.751 

2.837 

6.338 

3.918 

3.040 

40* 

14.790 

7.976 

2.873 

19.210 

7.236 

3.622 

9.949 

6.255 

3.179 

60* 

20.420 

10.700 

4.961 

19.990 

10.220 

4.702 

16.190 

9.149 

9.997 

•0* 

26.760 

JJ.250 

0.796 

24.240 

10.150 

5.605 

20.990 

10.760 

4.961 

Onat  and  Prager  [1]  gave  the  computation  formula  for  the  limit¬ 
ing  loading  capability  of  a  spherical  shell  without  opening  in  the 
presence  of  external  boundary  solid  support,  which  corresponds  to  the 
situation  in  this  paper: 

""2^(i-co»yj»Ci-Do[^(1~2K),<in,y,+^,(1~2K~4cMy,)cosyt  7) 

+  (3/C  — 1)  — ^^-(1—2/0' In  (1-2*)] 

where 


4J? 


Li  Kangxing^t  al.gave  an  even  more  convenient  formula  [2] 

w-sfl+je-j* )  (2-8) 

\  1—  cotp,  / 

In  order  to  compare  the  reliability  of  the  results  in  this  paper, 
a  special  case  within  the  range  of  applicable  conditions  is 

compared  to  (2-7)  and  (2-8).  The  numerical  results  are  shown  in  Table 
2-2  (equation  (2-2)  is  written  in  a  dimensionless  form  as  (2-7)). 


If  the  external  boundary  is  considered  as  the  simple  support  and 
the  internal  boundary  is  free,  then  by  using  an  identical  method,  we 
can  find  that  when  ,*•>*, 
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Figure  2-2 

sin  *<p,  (l-K'isia  *g>*—  sin  'y.  +  Ksin  V,+2(cosy*- cosy,)cos<p, 
3/f  (  cos  q>  i  —  cos  <p,)(  cosip,—  cosq>t+  sin 


when  **<P. 


cos'tp, _ tin  *y>~  «in  Vi _ 

*“  3  (coi^i  —  cosf>,)(cosf>,—  cosp,+  sin  Vi) 


(2-9) 


(2-10) 


where  9*»K  and  n  are  the  same  as  in  (2-^0,  (2-5)  and  (2-6) 


Reference  [3]  studied  the  opened  spherical  top  cover  under  the 
same  boundary  condition  by  using  the  "dual  moment  weak  interaction" 
yield  condition.  However,  its  expression  and  applicable  region  a  re 
very  complicated.  A  comparison  of  the  results  of  this  paper  to  those 
in  [31  is  shown  in  Figure  2-3- 


2.  Spherical  cover  top  with  a  center  opening  to  sustain  axially 
concentrated  load.  The  mechanical  model  of  the  structure  is  shown 
in  Figure  2-4.  When  reading  the  limiting  state,  let  us  assume  that 
the  distortion  speed  is 

Jfr~A(cosqsl  —  coi^j)  (2-11) 

where  is  the  position  of  the  plastic  region  which  is  to  be  deter¬ 
mined. 


Using  the  mobile  method  and  the  packed  shell  yield  condition, 
we  can  get 

2j*T  2K  +(1  — /Ocos1^,— 2cos$>,cos<p,  +~\~j^~c0>tP (2-1 


where 


ros  —  co*  <p, 


cos <p, -  cos*, - *in  <pt 


(2-12) 


(2-13) 


(2-14) 


(2-15) 


If  we  assume  that  the  structure  is  primarily  designed  to  sustain 
axially  concentrated  load  and  simultaneously  to  stand  the  uniformly 
distributed  internal  pressure  q,  the  mechanical  model  of  the  structure 
is  simplified  as  shown  in  Figure  2-5. 


If  we  assume  that  the  internal  pressure  does  not  change  the  dis¬ 
tortion  speed  (2-11)  under  the  limiting  condition,  we  can  obtain  the 
following  using  the  similar  method: 

.  r.  >  /«  « _  _  i  27C _ «  i 


where 


2«j ~ 2K  +  (1  —  K’)co5,f»i~2co«ylco«y,  +  cos'<pt  I 

cos  —  co*  <p, 

+  jjQ  (cosy t~  cos yt)*  +  ( cosy,  —  cos<pt)sin 
cosflp,  —  cos  9?, 

COS  9,—  COS  —  (cos 

Q- 

W  a.h 
a,Rh 

IK— (2  —  IK  +0)  cos  >,  —  20  sin  >cos< 

/"r  .  2(1-8*) 


(2-16) 


(2-17) 

(2-18) 

(2-19) 


(2-20) 
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TABLE  2-3 
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wawus* 
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1 

1  *  *<$ 

r, 

1  *  *<£ 

9 

MMS  (1)  © 

*  'U 
,  * 

* 

U  V 

1 

7«t«4D 

478 

6 

78 

6 

15 

r  V* 

15.4 

3.8 

2 

478 

8 

76 

0.6 

82.5 

33.3 

3.6 

3 

478 

8  - 

76 

8. 6 

32.0 

33.8 

4.8 

1 — specimen  no.;  2 — yield  stress  of  the  material  a  =28  kg/mm2; 
3,4,5 — mm;  6 — kg/mm2;  7 — limiting  axial  pressure  (fon);  8 — 
experimental;  9 — calculated;  10 — error;  11 — stainless  steel; 

12 — stainless  steel;  13 — stainless  steel; 


Equations  (2-12)  and  (2-l6)  agree  better  with  the  simulated 
tests.  The  comparison  of  the  experimental  results  to  the  theoretical 
computation  is  shown  in  Table  2-3. 


3.  The  triple  centered  top  cover  with  an  opening  in  the  center 
to  sustain  the  uniformly  distributed  internal  pressure.  The  triple 
centered  top  cover  is  made  of  part  of  a  spherical  shell  and  a  ring 
shell  which  are  smoothly  connected  together.  The  simplified  mechan¬ 
ical  model  is  shown  in  Figure  2-6. 


Let  us  assume  that  the  distortion  velocity  is: 
for  the  ring  shell  part 

W -A(l—  cin?) 

for  the  spherical  shell  part 

+  5(*in 
—  *in  <p) 

<{>  =  the  first  order  derivative  at  8. 

From  continuity,  we  get 


The  results  obtained  finally  are  tabulated  in  Table  2-4. 

where  r 

a,  —  arcsin  — 


(2-21) 

(2-22 


(2-23) 


26 


TABLE  2-4.  The  computation  formula  of  limiting  load  capability  of  a 
triple  centered  top  cover  with  a  large  opening  under  uniform  pressure 


TABLE  2-5.  Table  of  numerical  value  of  the  load  capability 
of  a  triple  centered  top  cover  with  rJD— 0 
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0.0600 

0.0740 

0.0730 

0.0903 

0.0001 

0.1 

0.20 

0.0807 

0.0403 

0.0017 

0.0004 

s  P/O. 
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\  7./0 

0.030 

0.034 

0.040 

0.040 

0.060 
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m 

HI 
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ft 

1 

0.16 

0.0006 

222 

0.0910 

0.0915 

0.1001 

0.1302 

1 

0.25 

0.0060 

0.0934 

0.0987 

0.1000 

0.1250 

0.1343 

0.1843 

0.1571 

1 

0.30 

0.0990 

0.0090 

0.1132 

0.1130 

0.1335 

0.1333 

0.1538 

0.1530 

BSS 

0.1073 

1 

0.35 

0.1053 

0.1055 

0.1106 

0.1197 

0.1410 

0.1411 

0.1024 

0.1024 

0.1760 

0.1770 

1 

0.40 

0.1112 

0.1110 

0.1201 

0 .1260 

0.1400 

0.1434 

0.1712 

0.1710 

0.1003 

0.1061 

0.0 

0.20 

9.0000 

0.0080 

0.1056 

0.1240 

0.1245 

0.1430 

0.1565 

0.1505 

2 — author;  3— reference  [4];  4 — this  paper;  5 — reference  [4j; 

6— this  paper;  7— reference  [4];  8— this  paper;  9— reference  [4]; 

10 — this  paper;  11 — author;  12 — reference  [4];  13 — this  paper; 

ill — reference  [4];  15 — this  paper;  16 — reference  [4];  17 — this  paper 

18 — reference  [4];  19 — this  paper;  20 — reference  [4];  21 — this  paper 
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Reference  [4]  completed  the  load  capability  computation  of  a 
triple  centered  top  cover  without  opening.  The  zero  center  opening 
angle  is  a  special  case  in  t^c  applicable  range  of  this  paper.  The 
comparison  of  this  work  and  reference  [4]  is  shown  in  Table  2-5. 

The  effects  of  various  opening  angles  on  the  load  capability  of 
the  triple  centered  top  cover  are  shown  in  Tables  2-6  and  2-7. 


TABLE  2-6.  The  effect  of  various  opening  angle  on  the  load 
capability  of  the  triple  centered  top  cover 
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TABLE  2-7.  Comparison  of  calculated  and  experimental 
results  of  load  capabilities  of  a  triple  centered  top  cover 
sustaining  uniform  internal  pressure 


r  i 

TT 

r, 

nr 

k 

~D~ 

R 

D 

P!°.  (f ) 

1 

0 . loss 

0.1000 

0.0002 

0.0770 

0.0200 

0.0200 

2 

0.1740 

0.1000 

0.0002 

0.9770 

0.0190 

0.0194 

4 — specimen  no.;  5 — material  a  =  *10  kg/mm2;  6 — load  capability, 
experimental  value;  7 — load  capability,  calculated  value;  8 — stain¬ 
less  steel;  9 — stainless  steel 

The  computational  formula  given  in  this  paper  was  proven  to 
satisfactorily  meet  the  engineering  approximation  requirement  through 
experimental  verification.  The  strain  plate  was  used  in  the  measure¬ 
ment  of  distortion.  The  theoretical  computation  can  be  simply  ob¬ 
tained  using  an  interpolation  method  from  Table  2-6  and  Figure  2-7- 
The  experimental  results  are  compared  to  the  computed  results  and  the 
comparison  is  shown  in  Table  2-7.  The  experimental  results  are  given 
in  Figures  2-8  and  2-9. 

4.  The  limit  analysis  of  the  severed  conical  top  cover  sus¬ 
taining  concentrated  load  and  the  minimum  weight  design. 

(1)  Equal  thickness  severed  conical  shell.  The  simplified 
mechanical  model  is  shown  in  Figure  2-10. 

The  limiting  load  capability  is  determined  by  the  following  for¬ 
mulas 

(* + r#)  coo  a+ oin ' a] «*■'«+  **+ r,(* + r,) ] 

L  ,  4  J  2  A_ri*  •  3  (2-2H) 

Jj-.'  + j(c«-  -iTj.),'- (Sr.c«,)  «  +  -«)  ♦  0]-«  o.jj) 


The  value  of  x  determines  the  position  of  the  plastic  zone  in  the 
limiting  state.  It  is  represented  by  the  distance  from  the  shell 
surface  to  the  rotational  axis.  From  experiment  and  analysis,  we 


Figure  2-9 


2 — strain;  3 — negative  value;  — negative  value;  5 — theoretical 
value;  6 — specimen  2;  7 — schematic  diagram  of  the  strain  position; 

9 — strain;  10 — specimen  1;  11 — theoretical  value;  12 — negative  value 
13 — schematic  diagram  of  the  strain  position 


Figure  2-11 

4 — linearly  thickening  severed  conical  shell;  5 — limiting  state 


know  that  the  destruction  of  the  severed  cone  caused  by  the  axially 
concentrated  load  only  occurs  in  the  local  region  in  the  small  end. 

If  the  plastic  zone  can  be  extended  toward  the  direction  away  from 
the  small  end,  the  loading  capability  can  be  improved.  Let  us  consi¬ 
der  such  a  conical  shell  thickness  variation  rate  (or  radially  ribbed 
structural  parameters).  The  limiting  state  corresponds  to  the  situa¬ 
tion  that  the  plastic  zone  occurs  in  the  entire  region  between  the 
two  end  planes  of  the  severed  cone.  The  material  utilization  rate  of 


this  structure  is  the  highest.  It  is  called  the  ideal  thickening 
severed  conical  shell. 

(2)  The  minimum  weight  design  of  a  severed  conical  shell  sus¬ 
taining  axial  load.  The  mechanical  model  of  the  linearly  thickening 
severed  conical  shell  is  shown  in  Figure  2-11. 


The  ideal  thickness  variation  rate  equation  is 

+ 30r*  8r5  W  +  [-^ir£r  cosa 

,cl2  *in‘a  .  .i  rth‘,  sin  ‘a  .1  q  ,nr,  ,  .  „ 


(2-26) 


To  satisfy  equation  (2-26),  the  load  capability  of  the  severed  con¬ 


ical  shell  is 


P-  *a.j[a(/?  +  r.)  +  ±b(R' +Rr,  +  r'i  )  ]cosa+(Af  -hi  ) 


(2-27) 


+fa.  tin  +y*9C*‘  +  r.(/?  +  r.)  ] 


in  (2-26) 


A—  1  —  - 


(2-28) 


The  mechanical  model  of  a  radially  ribbed  severed  conical  shell 
is  shown  in  Figure  2-12. 


The  ideal  reinforcing  rib  distribution  should  satisfy  the  follow¬ 
ing  parametric  equation: 

+£r-> 

(2-29) 

+2(»A+  - ?^pM3tf  +  rt)-*Ar;  +[aA*r.+»d// 

cma  9  y, 


TABLE  2-8 


2— specimen  no.;  3 — property  or  the  severed  cone;  4~materlal;  5— os  kg/mm’;  6,7,8,9,10, 
11,12— ran;  13 — q  kg/mm’;  14— experimental  value,  ton;  15— theoretloal  value,  ton;  16— equal 
thickness;  17— equal  thickness;  18— variable  thickness;  19— variable  thickness;  20— varlabl 

thickness;  21 — variable  thickness;  22 — equal  thickness;  23 — ribbed;  24 — ribbed;  25 ribbed; 

26— ribbed;  27— ribbed;  28— ribbed;  29—25 #  Bteel;  30—25#  steel;  31—25#  steel;  32— stain¬ 
less  steel;  33— stainless  steel;  34— stainless  steel;  35—25#  steel;  36—25#  steel;  37—25# 
steel;  38—25#  steel  39—25#  steel;  40—25#  steel;  41—25#  steel 


Figure  2-12 


To  satisfy  (2-29),  the  structural  load  capability  is 

P-«a,{[a(/?+r.)+  ^6(/?‘+7?rt  +  r;  )]co»a+ -*>"  *«  j  + 

♦ 

+*«#(//+*)  — +-7r-)+y9C/?‘  +  r.(^+r.)l 


(2-30) 


where  n  is  the  number  of  ribs,  H  is  the  height  of  the  rib, 
thickness  of  the  rib. 
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6  is  the 


The  reliability  of  the  above  analysis  can  be  verified  by  simula¬ 
tion  tests.  The  comparison  of  the  experiment  and  the  theoretical 
analysis  is  shown  in  Table  2-8. 


III.  THE  LOAD  CAPABILITY  OF  A  CENTRIFUGAL  TYPE  INJECTOR  DISK 


The  centrifugal  injector  disk  is  a  complicated  two  layer  circular 
plate  structure.  Its  double  layered  wall  divides  the  thrust  chamber 
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Figure  3-la 

2 — combustion  gas  hot  surface;  3 — bottom;  4 — nozzle;  5 — top; 

6 — combustion  gas  pressure  difference 

into  the  combustion  chamber  (see  I  and  II  in  Figure  3-1),  and  the 
oxidant  chamber  (II  or  I)  and  the  combustable  gas  chamber  (III).  The 
two  layer  plate  was  welded  together.  Generally  speaking,  the  pressures 
P1  and  P2  excited  on  I  and  II  are  approximately  equal  to  and  higher 
than  Pjjj-  Therefore,  the  injector  disk  of  the  entire  unit  only  has 
to  sustain  the  pressure  drop  Pj-Pjjj* 

The  mechanical  model  of  the  injector  disk  is  shown  in  Figure  4-16. 
In  the  computation,  the  following  basic  assumptions  are  adopted  [5]: 

i 

(1)  The  double  layer  circular  plates  were  connected  rigidly  by 
many  nozzles.  There  is  no  relative  distortion.  It  is  also  considered 
that  the  nozzles  are  uniformly  distributed  on  the  plate. 

(2)  The  bottom  of  the  injector  plate  is  working  at  high  temper¬ 
ature.  The  temperature  varies  nonlinearly  along  the  radial  flow  radius 
of  the  circular  plate.  The  temperature  is  lowest  near  the  edge  of  the 
circular  plate.  The  temperature  is  the  highest  at  the  center  of  the 
circular  plate.  Let  us  assume  that  the  temperature  variation  is 

W-r.-cr.-rjp-  (3-D 
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• 1 1 


where  Tq  is  the  temperature  at  the  cen¬ 
ter  of  the  plate  and  T.  is  the  temper- 

K  p 

ature  at  the  edge  of  the  plate.  p  =  ^; 
n  is  a  constant  greater  than  1. 
r  is  the  radial  flow  radius. 

Figure  3-2  The  limiting  structural  analysis 

assumed  that  thermal  stress  did  not 
affect  the  limiting  load  capability.  The  temperature  only  affects  the 
yield  stress  value  of  the  material. 

If  the  temperature  difference  between  the  center  and  edge  of  the 
circular  plate  is  not  too  large,  the  difference  in  the  value  of  the 
material  yield  stress  can  be  obtained  using  a  linear  difference  method. 
Therefore,  we  can  write  the  following  linear  relation: 

a.m-a.tl-KT)  (3-2) 

where  a — yield  stress  at  constant  temperature  (T=0) 
s 

K  — a  constant  which  varies  with  the  material. 

Consequently,  we  can  obtain 

<7»(r)  —  a,(e+rfp") 

where 

c-l-KT. 

(3)  The  yielding  bending  moment  of  the  double  layer  circular 
plate.  Under  the  yield  condition,  the  stress  distribution  along  the 
height  of  the  plate  is  as  shown  in  Figure  4-2.  Because 
a,(D-<r,(c+rfp'),  ot  =  constant. 


(3-3) 

(3-4) 
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(*l)  The  effect  of  the  nozzles.  A  series  of  densely  located 

nozzles  would  weaken  the  load  capability  of  the  circular  plate. 

Reference  [5]  suggested  that  this  effect  be  included  in  the  internal 

power  of  the  corresponding  velocity  field.  This  is  a  very  complicated 

method.  However,  [6]  recommended  the  use  of  a  functional  coefficient 

to  take  this  effect  into  account : 

I-1-- j?~  (3-6) 

F  — the  area  of  the  opening  of  the  nozzles  to  the  plate, 
o 

F  — the  total  area  of  the  plate. 


The  injector  disk  is  formed  by  two  layers  of  plates.  On  top  of 
the  double  layer  plate,  the  apertures  of  the  nozzle  are  not  equal. 
The  effect  coefficient  should  be  correspondingly  noted  as 


t  —  i  Ft  i 

*'  *  ~FT 

t  —  1  —  Fn 

*•  1  FT 


(3-7) 


Therefore,  the  opening  of  holes  on  the  nozzle  to  the  circular  plate 


is  considered: 


fr.CrW.Cc+rfp*)*, 

0,z*O'-£t 


On  the  aforementioned  basis,  if  the  edge  support  of  the  injector 
disk  is  considered  as  a  simple  support,  the  load  capability  sustain¬ 
ing  uniformly  distributed  pressure  is: 

P-5^-{l+^(c+4r)-«,^(c'+i^-+14r)}  (3-8) 

where  the  symbols  are  the  same  as  in  (3-5),  and  a  o i  .  it 

a*-  it 


If  the  edge  support  of  the  injector  disk  is  considered  as  a 
fixed  support,  then 

+u0y  (c*  /&)--**■  ( £&-)] 

where  is  determined  by  the  following  formula: 


(3-9) 


ApHm+"+Bpl'+Gpl+'-Dp:+Ep\lnp>  +  Kpl-L-0 


(3-10) 
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2(2ii  +  1) 

“  2(2n  +  l) 
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n  +  1 

li  »»a/9c(a/?c  —  2y)  ->- 1 
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IV.  THE  LOAD  CAPABILITY  OF  THE  COMBUSTION  CHAMBER 

The  load  capability  of  the  combustion  chamber  is  the  load  capa 

bility  of  the  combustion  chamber  itself.  The  weak  part  of  the  cham¬ 
ber  is  in  the  cylindrical  section.  It  is  the  representative  section 
for  checking  the  strength. 

The  cylindrical  section  has  densely  located  reinforced  ribs 
along  the  circumferential  direction.  The  inner  and  outer  walls  are 
connected  by  wavy  boards  longitudinally.  Obviously,  this  is  a  dual 
directional  jacketed  cylindrical  shell  with  reinforced  ribs.  It  can 
sustain  the  combined  effect  of  the  uniform  internal  pressure  and  the 
axial  force.  Assuming  that  the  cylindrical  shell  is  supported  rigidly 
on  both  sides,  let  us  also  assume  that: 

(1)  a  high  pressure  coolant  is  passed  between  the  Inner  and  outer 

walls.  It  is  considered  not  to  be  related  to  the  load  capability. 

(2)  The  inner  wall  temperature  of  the  combustion  chamber  is  the 
average  value  along  the  axial  and  wall  thickness  directions,  i.e.. 
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Figure  4-1 

1 — with  reinforced  ribs;  2 — connected  by  wavy  board;  3 — shape  of 
the  wavy  board 

the  inner  wall  yield  stress  is  a  constant.  Because  the  inner  wall 
is  thinner,  and  the  work  is  done  at  high  temperatures,  the  material 
property  is  low.  It  does  not  contribute  greatly  to  the  load  capability. 
Therefore,  this  type  of  simplification  will  not  bring  about  large 
effects . 

The  mechanical  model  of  the  simplified  cylindrical  section  of 
the  combustion  chamber  is  shown  in  Figure  4-1. 

Based  on  the  above  assumptions,  using  the  force  method  and  the 
simplified  yield  condition,  we  can  get  the  load  capability 


where 


1  +<*»  + 


4— jjMi£+4gg 
29J+c*(l+a.)! 


(4-1) 


are  the  relevant  dimensionless  parameters  with  respect  to  cir¬ 
cumferential  reinforced  ribs  and  axial  reinforced  ribs,  res¬ 
pectively; 

is  the  axial  force  coefficient,  which  reflects  the  cylindrical 
shell  characteristics  of  the  half  opening  of  the  combustion 
chamber 

is  the  structural  parameter  of  the  cylindrical  shell 
is  the  limiting  value  of  the  dimensionless  uniform  pressure. 


They  are  equal  to  the  following,  respectively: 


h1,h2  are  the  thickness  of  the  inner  and  outer  wall  respectively 
R  is  the  average  radius  of  the  cylindrical  section  of  the 
combustion  chamber 
R  is  the  thrust  of  the  engine 

P  is  the  combustible  gas  pressure  in  the  cylindrical  section 
of  the  combustion  chamber 

0^,02  are  the  material  yield  stresses  corresponding  to  the  inner 
and  outer  walls  at  average  temperature 
L  is  half  the  length  of  the  cylindrical  section  of  the  combustion 
chamber 

B  is  the  wave  distance  of  the  wavy  board 

A 

Bq  is  the  spacing  between  the  reinforced  ribs  in  the  circumferen¬ 
tial  direction 

o^q  is  the  material  yield  stress  of  the  reinforced  rib  in  the  cir¬ 
cumferential  direction 

F  is  the  cross-sectional  area  of  the  wavy  board  within  a  wave 

A 

length 

F0  is  the  cross-sectional  area  of  a  reinforced  rib  in  the  circum¬ 
ferential  direction 

is  the  yielding  variable  moment  on  a  unit  length  of  the  jacketed 
layer  formed  by  the  inner  and  outer  walls,  including  the  wavy 
board. 


H  is  the  height  of  the  wavy  board,  which  is  also  the  distance  between 
the  inner  and  outer  walls. 
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Figure  4-2.  The  load  capability  of  a  ribless  combustion 
chamber 

2 — experimental  curve;  3 — calcualted  value;  4 — hot  state  experiment 
and  calculated  curves;  5 — cold  state  experimental  and  calculated 
c  urve  s 


o^x  is  the  yield  stress  of  the  wavy  board  material 


+b)s 

\  un  w,  / 


where  6  is  the  thickness  of  the  wavy  board. 

Equation  (4-1)  agrees  well  with  the  experimental  results.  Figure 
4-2  shows  the  comparison  of  the  experimental  curve  and  the  theoretical 
computation  of  the  load  capability  of  a  ribless  combustion  chamber. 

V.  THE  LOAD  CAPABILITYOF  TUBE  BUNDLE  NOZZLES 

The  simplified  mechanical  model  of  the  tube  bundle  nozzles  is 
shown  in  Figure  5-1.  The  exerted  combustible  gas  pressure  is  assumed 
to  be 

P~P (5-1) 

Solving  the  set  of  equations: 

y.(f.)  +gt(.9t)P 
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Figure  5-1 

2 — tube  bundles — reinforced  rib;  3 — single  tube 

where  the  functions  are  as  follows: 

-rRNlivcwv-  »ing>)+M5 

JiW^AtCosv  +  RNKvcosq?—  sin  95)—  M%  — — ^ — — 

/.(»>) -  Ax  cos  93  +  At  s»n  q>  - -j-^T - a'  +  ia - 1 

|  /.(«>) -Atcotp+-g- 

,  x  *  R-H *in*>  MS  R(.Rv-Rnnvco$v-Hq>sinsp) 
gX<p) - - R  sin  <p-H  ' 

f  v  tytfc  .  R'H coa 9  _  ^[(Z-faQ  say+acwy] 


where  the  structural  parameters  and  material  characteristics  para¬ 
meters  are 


A‘-  **-*<»-*»w 

r(^,(<**in2p.+2oos2g3«)  ,  RH{  sin  <p,— acos<p.)  1. 

L  2(a*+4)  "+~ - 1+^ - 

_r  R‘(2+otg<pt)  ,  RH  1 • 

L  a'+4a  (1  +a*)cos^,  J  * 


A.-- 


1 


Rcos<p, 


a  _f  i?*(2sui2g»,+2co4  29>.)  _l  RH(sin q>,— acos^p,)  ~\ 

4  L  2(a*+4>  +  r- /  * 


VI.  THE  LOCAL  STRENGTH  OF  THE  COMBUSTION  CHAMBER 

The  so-called  local  strength  is  the  uniform  load  resistance  of 
the  cooling  conduit  of  the  combustion  chamber  and  the  connection 
strength  of  the  wavy  board.  Their  characteristic  is  that  the  struc¬ 
tural  dimension  is  so  small  that  it  is  extremely  difficult  to  reason¬ 
ably  simplify  the  computational  model.  The  computational  formula 
given  here,  after  experimental  verification,  has  sufficient  accuracy 
for  the  specific  problem  of  the  engine. 

(1)  The  connection  strength  of  the  wavy  board.  The  simplified 
mechanical  model  is  shown  in  Figure  6-1.  Let  us  take  an  element  in 
the  combustion  chamber  for  analysis. 


Figue  6-1 

2 — Inner  wall;  3 — wavy  board;  ft — outer  wall 

Let:  a  Is  the  inner  wall  yield  stress  of  the  combustion  chamber 

aB  is  the  strength  limit  of  the  wavy  board  material 

o^drill  represents  the  strength  limit  of  the  drill  weld 
6  represents  the  thickness  of  the  wavy  board 
n  is  the  thinning  coefficient  of  the  wavy  board  after  the 
impulse  pressure  treatment 

<j>0  is  the  angle  between  the  side  of  the  wavy  board  and  the  cir¬ 
cumferential  direction  of  the  combustion  chamber 
Bx  is  the  wave  distance  of  the  wavy  board 
b  is  the  width  of  the  drill  weld 

h  is  the  inner  wall  thickness  of  the  combustion  chamber 
PMp  is  the  pressure  distribution  of  the  cooling  conduit. 


The  pressure  limit  inside  the  cooling  conduit  can  be  determined 
by  the  strength  limit  of  the  drill  weld  seam  as: 


drill 

at  o*  iny« 


B,  »in  97, — 2i 


(6-1) 


The  limiting  pressure  of  the  cooling  conduit  can  also  be  deter¬ 
mined  by  the  strength  limit  of  the  wavy  board  itself  as 


*  wav6  . 

p  _  Ztr#  'O’  stop,  _ 
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(2)  The  load  capability  of  the  cooling  conduit 

D  _  16  a.k* 


(6-2) 


(6-3) 
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The  drill  weld  strength  is  a  statically  determined  problem. 
Equations  (6-1)  and  (6-2)  are  the  destructive  strength,  which  is 
used  to  replace  the  load  capability. 

VII.  THE  RING  SHAPED  FLUID  COLLECTING  RING 

The  schematic  diagram  of  the  structure  is  shown  in  Figure  7-1. 

The  uniform  pressure  load  capability  [7]  is: 

„  \. - ! - 


1  +~na 

the  upper  limit  solution - f -  (7_1) 

l+_  na 

Fieure  7—1  P,  tr»  1  •4*o 

the  lower  limit  solution  down--^-  i 

1+—  *a 

(7-2) 

where 
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